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Physical spin configurations corresponding to topological excitations , expected to be 
present in the XY limit of a quantum spin | Heisenberg anti-ferromagnet, are probed on a 
two dimensional square lattice . Quantum vortices (anti-vortices) are constructed in terms 
of coherent staggered spin field components, as limiting case of meronic ( anti-meronic ) 
configurations . The crucial role of the associated Wess-Zumino-like ( WZ-like ) term is 
highlighted in our procedure . The time evolution equation of coherent spin fields used 
in this analysis is obtained by applying variational principle on the quantum Euclidean 
action corresponding to the Heisenberg anti-ferromagnet on lattice . It is shown that the 
WZ-like term can distinguish between vortices and anti-vortices only in a charge sector 
with odd topological charges. Our formalism is distinctly different from the conventional 
approach for the construction of quantum vortices ( anti-vortices ) . 
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1. Introduction 

One of the most convenient approaches to investigate the origin of topological exci- 
tations in quantum spin systems is through the coherent state formulation [1-5]. The 
geometrical phase ( Berry phase ) which arises in the partition function of the quantum 
spin models in this approach , plays the central role in deciding the topological sectors 
of various spin field configurations . For instance, in the case of the Heisenberg antiferro- 
magnetic chain in the long wave length limit, one has a two-dimensional ( one spatial plus 
one pseudo-temporal) nonlinear sigma model . This model contains the geometrical phase 
which acts as the topological WZ term [1] . For lattices with spatial dimensions of two or 
more, such type of term vanishes in the long wavelength limit [4]. Interestingly enough, 
this term is found to survive in the medium wave length limit and finds its signatures in 
various numerical and analytical calculations, as well as in neutron scattering experiments 
[6-9]. It may be recalled that in the long wavelength limit, one retains terms up to first 
order in spatial derivatives in the expansion of the spin field within the WZ-term. In this 
limit, the WZ-term survives in the case of a spin chain but vanishes for a square lattice 
[2,4]. Therefore in the latter case, it is necessary to retain terms up to 2nd order in spatial 
derivatives to extract a non-vanishing contribution. This requires a total length scale of 4 
lattice spacings (viz., 2a in the positive direction and 2a in the negative direction, a being 
the lattice parameter ) along both the axes. This translates in the q space as a regime 
around (l^l? |^|) • Thus we are in the vicinity of the mid-zone region in the q-space . 
Calculations based on this region are referred to as "medium wavelength approximation". 
Hence it is extremely important to perform calculations explicitly on the lattice in order 
to get a finite contribution to the WZ-term in the lowest order . The geometrical phase 
in its discrete form is the lattice analogue of the Wess-Zumino (equivalently Berry phase) 
term. Henceforth we refer to this discrete expression as "WZ-like" term . 
One of the major motivations for our present study has been theoretical determination 
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of the dynamical structure factor 5'(q, u) for low dimensional quantum antiferromagnets, 
particularly in view of the inelastic neutron scattering experiments on the cuprates [6] 
. In particular, the occurrence of 'central peak' in the experimentally observed ^(q, u) 
in the medium wavelength regime, strongly indicates possible existence and dynamics of 
topological excitations [6,7]. The 'central peak' refers to the peak occuring at a; = 0, 
in the plot of S'(q, u) versus u in the constant g-scan. Generally, this is an important 
signature for the spin dynamics driven also by the translational motion of the topologi- 
cal excitations and defects [3, 7]. Furthermore , the spin-spin correlation length ^(T) as 
calculated from the experimentally extracted C(q) ( the static q-space spin-spin correla- 
tion function ) for spin ^ two-dimensional quantum Heisenberg anti-ferromagnet (QHAF), 
shows marked departure from the renormalized spin wave theory, in the cuprates [9,10] . 
These observations combined with the theoretical analysis described above, suggest that 
by retaining the discrete lattice structure for low dimensional quantum antiferromagnets 
it should be possible to get a nonvanishing contribution of the Berry Phase or the WZ-like 
term [3,8,9 ] . 

Earlier detailed analysis was done by us for a general quantum anisotropic Heisenberg 
spin model using coherent state formalism in a quasi-continuum limit [3,8]. Now we deal 
with the extreme quantum case of a strongly anisotropic i.e., XY-limit of a spin | QHAF 
explicitly on a lattice. 

It is worthwhile to refer to another approach to study topological excitations [9] . In that 
approach the technique of Schwinger Boson Mean Field Theory (SBMFT) was applied on 
a quantum Heisenberg anti-ferromagnetic spin system at an intermediate length scale , 
assuming that the corresponding field theoretic action does not contain any Berry phase/ 
Wess Zumino ( topological like )term , in contrast to the real situation. Therefore the 
topological excitations are introduced heuristically there [9] . 

The main aim of this paper is to verify the presence of the WZ-like term and examine 
the properties of the physical spin configurations . We demonstrate that the WZ-like 
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term really identifies a large class of the topological excitations and clearly differentiates 
between quantum vortices and antivortices with different charges . 



2. Mathematical Formulation 
2.1 Action and the WZ-like term 

The quantum Euclidean action Se for the coherent spin fields n(r,t) can be written as 
[2] 

Se = -isgY,Swz[m{r, t)] + f dm{rv) (1) 

r Jo 

where 's' is the magnitude of the spin (s = ^ in the present case ) , r is the position vector 
of a lattice site and 

(n|S|n) = sn 
?/(n) = (n|H(S)|n) 
|n) = nin(r,t)) 

r 

(2) 

where denotes the direct product of all coherent spin states over the spatial lattice . 
The quantity 'H(S) is the spin Hamiltonian on the lattice and 5H/z[m(i"5 1)] is the WZ-like 
term on a single lattice site. 

'H(S) in our calculation is given by equation (5) in spin s representation and 'H(n) the 
corresponding Hamiltonian in the coherent spin fields n(r,t) ( see equation (6) ). The 
variable t denotes the pseudo-time (Euclidean time) appropriate to the above fields and 
has the dimension of inverse temperature . The parameter (3 stands for ^ as usual; T 
being the real thermodynamic temperature of the spin system. The WZ-like term Swz 
corresponding to a single spatial lattice point r at a fixed time 't' is given as follows [2]: 

5vi/z[m(r,^)] = dt (irm(r, t, r) ■ 9tm(r, t, r) A 9rm(r, t, r) (3) 

JO JO 
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with m(r,t,0) = n(r,t), m(r,t, 1) = no(r), and m(r, 0, r) = m(r,/3,r), t e [0,(3], 
T G [0, 1] . 

The expression in equation(3) is the area of the cap bounded by the trajectory F parametrized 
by n(r,t) [={ni{r,t),n2{r,t),n3{r,t))] on the sphere: 

n(r,t)-n(r,t) = l (4) 

Furthermore the fields m(r, t, r) are the fields in the higher dimensional (t, r)-space and 
the boundary values n(r, t) are the coherent spin fields . The field no(r) is the fixed point 
(0,0,1) on the above sphere and the state vector |n(r,t)) appearing on the right hand 
side of equation (2), is the spin coherent state at a single lattice point r [1-3] . 
For the two dimensional lattice , we express n(r,t) = n{ia,ja,t) as n{ia,ja) for breavity 



2.2 The Spin Hamiltonian 

The spin Hamiltonian corresponding to anisotropic Heisenberg spin system of XX Z type 
with antiferromagnetic coupling is given by 

H{S) = g Y: S(r) ■ S(r/) + gX, ^ S.{r)S.{rf) (5) 

(r,r/) (r.r') 

with g > and < < 1 , r,r/ running over the lattice, and (r, r/) signifies nearest 
neighbours and S = (S,S'2). Here S = {Sx,Sy) is the projection of the operator S onto 
the XY-plane . 

It follows from equations (2) and (5) that the spin Hamiltonian in terms of coherent spin 
fields is given by 

Hin) = g s"^ h{ia,ja)-n{i'a,j'a)+gs^Xz ^ n^{ia, ja)n2.{ra, so) (6) 

Here n(m,ja) = {ni{ia, ja),n2{ia, ja)) is the projection of the coherent spin field 
n onto the XY-plane . Mathematically , n{ia,ja) = (n|S(ia, ja)|n) where the state 
vector |n) is given by equation (2) . 
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From equations (1), (2) and (6) the quantum action (Euclidean) for the two dimensional 
anisotropic spin system with antiferromagnetic coupling is given by [2] 

Se = —is^^Swz[™-{ici, ja)] + / dt[ g s'^ ^ n{ia, ja) ■ n{i'a, j'a) 

+g\zS^ n-i{ia,ja)n^{ia,i'a)] (7) 

((i,i),(*',i')> 

where the constraint given by equation (4) has to be satisfied at each lattice point . 
As our system is a quantum antiferromagnet , we assume the spin configurations to have 
correlations exhibiting a bipartite symmetry , when the temperature is not too high [8]. 
We therefore stagger the configuration as: 

n(m,ja) — > (-l)*+^n(za, ja) (8) 

Thus we obtain from equations (4) , (7) and (8) the following total quantum Euclidean 
action for the anisotropic Heisenberg antiferromagnet 

ST'' = -^sY: i-^y^' Swz[ni{t.a,ja)] + dt g I i-iy''''"'^''n{ia,ja) ■ fi( 

- dt Yl^^^Kj (^'2' i«) - 1] 

Hence for the nearest neighbour interaction the above action is given by: 

ST'' = -^sY: i-iy^' SwzlMhj)] + fdt{ - g .T[n(z, j) ■ n(z + 1, j) 

+n(z,j) ■ n{i,j + 1)] - g X^s'^YMh Jhsii + 1, j) + n^ii, j)n3{i, j + 1)] 

(id) 

-f^dtY^'KA^'ihj)-!]} (10) 

The last term on the right hand side of the above equation is the lattice version of the 
term /(i^x/o^A(x, t)(n^(x, t) — 1) where A(x, t) is an auxiliary field playing the role of a 
multiplier . The equations of motion which follow from the minimisation of the above 
action are given by equations (A.l) in the Appendix . The purpose of this exercise is 
to obtain an expression for S^"^^ given by equation (11) below , in terms of individual 
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contributions from lattice sites. The quantity S^^^ is topological-like in the sense that 
the WZ term survives in the medium wave length limit and it indicates topological ex- 
citations [1, 8-10] . One can derive the expression for S^^^ by using equations (11)-(13) 
below ( for example , as we have done in the case of 1-vortex ;resulting in equation (A. 5) 
in the Appendix ) . Then we substitute for 9tn from the equations of motion [ A.l ) into 
the expression for S^^^ ( in equation ( A. 4 ) ) . 

The first term on the right hand side of the above Eqn. (10) is the lyZ-like term on the 
lattice which is formally written as: 

5^T = E(-ir^^'5M/z[m(^,j)] (11) 

It is possible to evaluate only the 'difference' of 5vi/z[ni(za, ja)] terms from two neigh- 
bouring lattice sites, in terms of the coherent spin fields n(za, ja). Thus to extract the 
topological-like contribution from S^^^ , we use the following expressions for the above 
'difference' [2] : 

5a;5vKz[m(r)] = Swz[m{ia,ja)] - Swz[m{{i - l)a, ja)] 

3 

dt[(5,n- (nAatn)](r) 
5ySwz[T^{r)] = Swz[m{ia,ja)] - Swz[m{ia, {j - l)a)] 

= rdt[6yn-{nAdtn)]{r) (12) 

J 

where 

6xii{ia,ja) = n{ia, ja) — n{{i — l)a, ja) 

Syn{ia,ja) = n{ia, ja) — n{ia, {j — l)a) (13) 
with r = {ia,ja) , and i,j = 1, 2, , 2A^ . 

From equations (11) - (13) we can write S^^^ as given by equation (A. 4) in the Appendix 







3. Calculations and Results 
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We now construct vortices ( anti-vortices ) in the flattened meron conflguration limit with 
staggered spin flelds . It is important to point out at this stage that it is necessary to 
retain both and as inflnitesimal but nonvanishing quantities , in order to generate 
a topological character of the above configurations ( see equation (7) and Ref. [3] ). Af- 
ter the construction of an elementary vortex ( anti-vortex ) plaquette with coherent spin 
field components ni{ia,ja), n2{ia,ja) and n^iia.ja) we present calculations for S^^^ 
corresponding to vortices (anti- vortices). We analyse all the higher charged vortices ( 
anti-vortices ) by decomposing them in terms of elementary plaquettes of unit charge . 
Our findings are : 

I) Odd-charged configurations can be described consistently in our scheme . 

II) Even-charged configurations however , behave in an anomolous way , which is ex- 
plained in the Appendix . 

III) Simple calculations based on the Hamiltonian with staggered spin fields [ see equa- 
tions (10) and (A. 6)] show that the excitation energy of a 1-vortex, as measured from 
the ground state, is approximately +12gs^ . It may be noted that the ground state spin 
configuration is obtained by taking the coherent spin vectors (n(za,ja)) aligned in the 
same direction . These topological excitations can indeed be produced even at very low 
temperature by means of quantum fiuctuations , thereby maintaining a steady vortex- 
anti-vortex pair density even at temperature close to zero [14] . 

3.1a Analysis of 1-vortex 

We assign the coordinates ((z — l)a, (j — l)a), (m, (j — l)a) , {ia,ja) and ((z — l)a,ja) 
to the vertices 1, 2, 3 and 4 respectively, of an elementary 1-vortex [Fig.l] plaquette and 
construct the quantum vortex with staggered spin fields. 

As we are interested in the extreme XY- anisotropic limit, we assume n^{i'a,j'a) = 
sin e{i'a, j'a) at each lattice point {i'a,j'a) , where e{i'a,j'a) is a very small positive quan- 
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Fig. lb 



Figure 1: (a) 1-vortex , (b) 1-anti- vortex 

tity dependent on lattice site {i'a,j'a). Then it follows from equation (4) that within a 
vortex (or anti-vortex) the following solutions emerge : 



At any vertex , rii = ±(1 — S) ; n2 = iy [cos^e(/a,ma) — (1 — 5) ] and vice versa . 



For illustration [ Fig.-l ] we have a quantum vortex of charge +1 in which the horizontal 
arrow — ^ at a vertex represents rii with value 1 — 6 and the vertical arrow f implies n2 
with value 1 — 6. Further, in this figure the horizontal arrow at a vertex denotes rii 
having value —{1 — 6) and the vertical arrow l represents n2 with value —{1 — 6). For 
vortex with stable spin configuration the horizontal or vertical arrow should have a steady 
magnitude of ±(1 — 5) . This implies that the quantity 6 must be independent of time ( 
described in the Appendix ). It may be remarked that the vortex configuration (shown 
in Fig. la) has a topological charge +1 , as the spin rotates through an angle +27r in 
traversing the boundary once in the anti-clockwise sense. 

In this connection let us point out that usually in a two dimensional vortex corresponding 
to spin i quantum spin model, the states | — )■) , | ) are taken to be the eigenstates of 
with eigenvalues +| , — | respectively. Similarly the states | t) ) I i) ^^^^ the eigenstates 
of Sy with eigenvalues +| , — | respectively [11] . However here we are not following that 
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scheme . We make use of the coherent states given by |n) = cos f ||) + (e~''^)sin 2 ^'^^ 
evaluate the expectation values of and Sy (referred earlier as rii and n2 respectively ).In 
our picture the horizontal and vertical arrows represent these expectation values rii and 
n2 respectively . To be more precise our configurations are in fact 'flattened merons(anti- 
merons)', mimicking a vortex (or anti- vortex) in the limiting case — [12, 13] . 

This prescription will be followed for constructing the vortices {anti — vortices) 
of higher charge values as well . 

3.1b Contribution of S^^^ 

Let us now look at the symmetry properties of the contribution of S^^^ to the vortex 
plaquette, which we have denoted by [Swz\{i_.uortex) ■ shown that (see Appendix) 

the contriburion of iS^f ^ to a vortex plaquette can be decomposed into two parts such 
that the first part remains invariant when going over from a 1- vortex configuration to the 
corresponding anti-vortex configuration ; whereas the second part changes sign under this 
operation. We denote the first part by A and the second part by B. This transformation 
from vortex to anti- vortex is implemented by changing 722(2) and n2(4) in Fig.l to -^2(2) 
and -^2(4) respectively, in the case of a 1-vortex. Algebraically this means that A contains 
terms which are quadratic or of even degree in ^2(2) and n2(4) while B contains terms 
that are linear or of odd degree in n2(2) and ^2(4). Using equations (11)-(13) we cast 
Sy/'z^ in the form A + B as explained in the Appendix . 

We adopt the following algorithm for the construction 
of higher vortices (anti- vortices) :- 

For the construction of higher vortices we calculate the contribution of S^^^ given by 
equations (11)-(13) on a plaquette by algebraically adding the contributions of S^^^ on 
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each of the individual elementary plaquettes (subvortices) with a weightage factor of | to 
the common bonds shared between the pairs of adjacent subvortices. We are interested in 
those field configurations for which the contributions of S^^^ on the common bonds cancel 
each other (see Appendix) and only the contribution on the peripheral boundary remains. 



3.2 Analysis of 2-vortex 
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Figure 2: 2-vortex 

For a typical 2-vortex we refer to Fig.2 . We denote the subvortices by a , b , c and 
d, each carrying topological charge -|-1 [ Fig.2 ] . However the spin at the central lattice 
point of the vortex ( the point O in Fig.2 ) becomes non-unique , as is clear from the 
construction of the subvortices . Thus the central point situated on a lattice site turns 
out to be a 'defect' or a singular point and cancellation of the contribution of S^^^ along 
the common bonds does not lead to a consistent spin field configuration . In other words 
the contributions of S^^^ on the boundary of the 2-vortex is not well defined and the 
construction of 2-vortex in this scheme becomes problematic . The scenario persists in all 
the vortices (anti-vortices) possessing even valued topological charges , as can be read 
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out from the spin field configuration in the case of 4-vortex [F2(yf.4] (See Appendix ) 



3.3 Analysis of 3-vortex 

For 3-vortex we have a consistent spin field configuration [ Fig. 3] . We have now an 
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Figure 3: 3-vortex 



elementary anti-vortex plaquette at the central region with well defined staggered spin 
field configurations and the contributions of S^^^ along the common bonds cancel each 
other, giving rise to a consistent spin field configurations . It may be pointed out that 
in contrast to the case with even- valued charge , we now have a subvortex with opposite 
charge occupying the central region. 



3.4 Analysis of 4-vortex 



Similar to the situation corresponding to the 2-vortex , here again the spin at the 
centre becomes ambiguous, owing to the vanishing of the effective horizontal component 
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Figure 4: 4- vortex 

[ Fig. 4] . Thus the central point turns out to be a defect or a singular point . As a 
result the contributions of S^^^ on the common bonds do not cancel each other . 
The ambiguity of the suitable spin field configuration for 4- vortex is explained in the last 
paragraph of the appendix. 



4. Conclusion and Discussion 

i) Our calculations and analysis with 2D spin | anti-ferromagnetic quantum XY model on 
discrete lattice clearly brings out distingushing features between the even and odd charge 
sectors. It may be recalled that the presence of a topological term in a field theoretic 
action indicates that the field configurations can be classified into different sectors with 
various topological charges . 



ii) Our work has firmly established the role of WZ-like term , derived from a micro- 
scopic theory , as a topological charge indicator for the odd-charged excitations . Our 
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scheme is implemented in such a way that only a net contribution along the outer bound- 
ary of a vortex plaquette remains nonvanishing . As a result , the WZ-like term becomes 
an implicit function of the charge . 

iii) Different type of quantum mechanical operators were used to determine topological 
charges of excitations in spin models [14] . In our present work on QHAF we used a field 
theory based approach to identify such excitations. 

Our future plan includes the generalization of our approach to the case of finite to 
achieve a physical realization of meronic type of excitations in the spin models. As we 
increase the magnitude of Xz from infinitesimal to a finite value , the vortices/flattened 
merons would gradually go over to complete 3-dimensional merons . 
Furthermore , we aim to evaluate the static and dynamic spin correlations for two di- 
mensional spin I anisotropic quantum Heisenberg anti-ferromagnet at any temperature 
by taking into account the contributions from these topological excitations as well as 
spin waves. These calculations will be very much useful to critically examine the novel 
concept of "quantum Kosterlitz - Thouless transition" corresponding to quantum anti- 
ferromagnets [ 3, 15, 16, 17 ]. Above the Kosterlitz-Thouless transition temperature where 
the vortex-antivortex pairs unbind , the dynamical structure factor is expected to exhibit 
a central peak [6]. 

Our picture has an interesting analogy with the work presented in [18]. However there 
are some finer differences. In our case we only have electrically neutral topological exci- 
tations which occur at zero temperature as well as at finite temperature. On the other 
hand, in [18] there are both electrically charged and electrically neutral type of topologi- 
cal excitations. Whereas the charged ones appear at zero temperature, the neutral ones 
appear only at finite temperature . 
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To conclude, our study of topological spin excitations on the 2D-lattice will un- 
doubtedly play an important role in the analysis of thermodynamics of low dimensional 
ant i-f erromagnet s . 

Acknowledgement : One of the authors (SKP) would like to thank M.G.Mustafa, Ra- 
jarshi Roy and Purnendu Chakraborti for their valuable help in the preparation of the 
figures. 
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Appendix A 

Applying variational principle to the action given by equation (10), we obtain the 
equation of motion for an XY-anisotropic Heisenberg antiferromagnet explicitly on the 
lattice , in the following form: 

(-1)*"^^ dtni{ia,ja) = -i g s [A^risiVs - n3N2]{ia, ja) 
(-1)*"^^ dtn2{ia,ja) = -i g s [nsNi - X^niN3]{ia, ja) 

(-1)*"^^ dtn^{ia,ja) = -i g s [niN2 - n2Ni]{ia,ja) (A.l) 
In the above derivation we have used the following variation of the WZ-term given by [2] 
(55^1^[m] = J2i,j {-^y^^ 5Swz[m{^a,ja)] 

= (-1)'+' /o^ dt 6n{za,ja) ■ (n A dtn){za,ja) (A.2) 

In Eqn.(A.l) Ni{ia, ja), N2{ia, ja), N^lia, ja) are the components of the vector N at the 
lattice point {ia,ja) . The vector N is given as: 

N{ia,ja) = n(za, (j — l)a) + n((z — l)a,ja) + n((z + l)a,ja) + n(za, (j + l)a) (A. 3) 
From equations (11), (12) and (13) we can write S^^^ in the following form : 

2Spi' = 2E,, (-l)^"'^' SwzHta,ja)] 

= E,, i-iy-'Io^ dt {[n(z, j) - n(z, J - 1)] ■ (n A 9,n)(z, j) 

+ [n{i,j) -n(z - l,j)] ■ {nAdtn){i,j) - [n(z,j - 1) - n(z - l,j - 1)] ■ {nAdtn){i,j - 1) 

-[n(^ - 1, j) - n{t - 1, J - 1)] • (n A dtn){z - 1, j)} (A.4) 
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Thus for the 1- vortex ( anti- vortex ) plaquette having vertices 1,2,3 and 4 (see Fig.l 
) with coordinates ((i — l)a, (j — l)a), (m, (j — l)a), {ia,ja) and {{i — l)a,ja) respectively 
the contribution of S^^^ to the 1-vortex is given by 

[Swl%-.,orte.) = -Io'dt[n{{z - l)a,ja) ■ (n A {-iy^^d,n){za,ja) 

+n(m, (j-l)a)-(nA(-l)*^^9tn)(m, ja) + n((i-l)a, (j-l)a)-(nA(-l)*"^^9tn)(m, (j-l)a) 

+n((2 - l)a, (j - l)a) ■ (n A dtn){{z - l)a,ja)]} (A.5) 

We evaluate the right hand side of Eqn.{AA) by substituting for dfH from Eqns.{A.l) and {A.3) 
. To keep the calculations simple but consistent , we retain the intra-plaquette contribu- 
tions by imposing a so called "local periodic boundary condition" (local PBC) as applied 
to the site closest to the vertices belonging to the plaquette under consideration . 
For example, in the case of 1-vortex [see Fig. la] the local PBC implies: 

n(m, ja) = n((i — 2)a,ja) = n(m, (j — 2)a) 

n((« - l)aja) = n{{i - l)a, (j - 2)a) = n((z + l)a,ja) 

n(m, (j - l)a) = n{ia, (j + l)a) = n((z - 2)a, (j - l)a) 

n((z - l)a, (j - l)a) = n((2 + l)a, (j - l)a) = n((^ - l)a, (j + l)a) (A.6) 

Besides we make use of the following conditions satisfied at different vertices of the pla- 
quette for all time, as explained in the section 3.1 [see Fig. la] : 

ni((i — l)a, (j — l)a) = —ni{ia, ja) = 1 — 5 

and 

n2(m, (j — l)a) = — n2((i — l)a,ja) = 1 — 6 (A. 7) 

This means that the following equations must hold : 
dtni{{i - l)a, (j - l)a) = -dfUi^iaJa) = dt{l - 5) 
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dtn2{ia, (j - l)a) 



dtn2{{i - l)aja) = dt{l - S) 



(A.8) 



Making use of Eqns 



(A.3) and (A. 6) in (A.l) we have 



dtni{{i 



i g s \^n2N^{{i - l)a, (j - l)a) 



i g s X:,n2N3{ia,ja) 




dtn2{ia, (j 



l)a) = i g s X^niN-s^ia, {j - l)a) 



{-ly ^^^dtn2{ii - l)a,ja) = i g s X;,niN3{{i - l)a,ja) 



(A.9) 



The right hand sides of Eqs. (A.9) vanish in the flattened meron configuration hmit where 



ration. Consequently from Eqns. (A.8) it follows that the quantity 5 does not vary with 
time . Note that 6 is independent of the position of the vertex points of the plaquette 
[See Fig.l . clj , clS explained in section 3.1 

It is important and interesting to point out that had we not used the local PBC (A. 6), we 
would not have had "time independent 5" . In other words, the spin field configurations 
forming the vortices become non-stationary in absence of the local PBC . For example 
using the 1st equation under (A.l) and the equation (A.3) we have : 



(-1)' ^'^^ ^dtni{{i - l)a, (j - l)a) = -i g s [\zn2Ns{{i - l)a, {j - l)a) - n2i{i - l)a, (j - 
2)a)-n2((^-2)a, (j - l)a)] 



The right hand side ( rhs ) of the above equation does not vanish even in the flattened 
meron configuration limit. This is unlike in the case of the 1st equation under (A.9) where 
the parameter 6 came out to be time independent. However application of the local PBC 
(A. 6) restores the rhs (as referred to earlier) to that of the 1st equation under (A.9). 

Now using the Eqns. (A.l), (A.3), (A. 5), (A. 6) and (A. 7) we obtain Swz as follows : 



and N3 acquires very small value , since e{ia,ja) is very small in this configu- 



{I— vortex) 



[ni(2) + ni(4)][(n ■ N)(3)ni(3) - iVi(3)] 
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+ [n3(2) + n3(4)][(n ■ N)(3)n3(3) - N,{3)] + n,{l)[{n ■ N)(4)ni(4) + (n ■ N)(2)ni(2)] 
+n2(l)[(n ■ N)(4)n2(4) + (n ■ N)(2)ri2(2) - 2iV2(2)] + ns{l)[{n ■ N)(4)n3(4) 
+ (n-N)(2K(2)-2iV3(2)] 

= A + B (A.IO) 

where in the flattened meron configuration hmit A and B are given by: 
A = j^dt {2(ni(2) + ni(4))'[ni2(3) - 1] + 2 ni(3) ^3(3) [ni(2) + ni(4)][n3(2) + ^3(4)] 

+4 ^2(1) [V(2)-1]K(3) +722(1)]} 

B = Jl^dt {+2 ^2(2) ni(l) [n,{2) - m{A)] [n^iS) + ^2(1)] 

+ 2 ^2(2) 713(1) [^3(2) - ^3(4)] [712(3) + 712(1)]} (A.ll) 

In the equations (A.IO) and (A.ll) above we have labelled the vertices of the 1-vortex 
[Fig. la] by 1, 2, 3 and 4 in the anti-clockwise sense. The quantity 77j(l) denotes the i-th 
component (i = 1, 2, 3) of the coherent spin vector n at the vertex 1 and similarly 7ij(2), 
77 j (3) and 7ij(4) stand for the i-th components at the vertices 2, 3 and 4 respectively . 
It is interesting to note that A remains invariant if we go from vortex to antivortex by 
changing 72,2(2) and 72,2(4) in Fig.l to —722(2) and —772(4) respectively whereas B goes 
over to —B . Thus [<Sw'z^](i-vortex) t^-kes the form A — B ioi antivortex [Figl.b]. 

To analyse the case of 2- vortex [Fig. 2] we assign the co-ordinates to the vertices 1, 
2, 3, 4, 5, 6, 7, 8 as ((7 - l)a, (j - l)a) , (m, (j - l)a), ((7 + l)a, (j - l)a), ((7 + l)aja), 
((7 + l)a, (j + l)a), (m, (j + l)a), ((7 — l)a, (j + l)a), ((7 — l)a,ja) respectively . The 
centre O has co-ordinates (ia,ja). Using Eqns. (A.l) , (A. 3) as well as the local PBC 
, as appropriate to the case of 2-vortex , the equations of motion of the spins 712 at the 
vertices 2, 4, 6, 8 reduce to the following , in the flattened meron conflguration limit ( 
A, ^0): 

20 



(1 - Sf = cos\{2) = 005^6(4) = cosh{6) = cos\{8) (A.12) 

It is therefore obvious from the above equations that the z-component of the spin n viz., 
77,3 becomes position independent. 

In Fig.3 we have denoted the elementary vortex plaquettes by a, b, c, d, e, f , g, h. Note 
that there is an elementary anti- vortex plaquette in the central region with vertices A, 
B, C and D . We write formally the expression for S^^^ f^^_^^^^^^^ : 

\ostagg^ _ 
V-^WZ i(3-vortex) — 

J^dt{n{la)inAdtn){2a) + ln{2a)inAdtn){3a) + ln{4a)inAdtn){3a)+n{la)inAdtn){4a) 
+n(16) ■ (n A dtu) (26) + n(26) ■ (n A dtu) (36) + |n(46) ■ (n A dtu) (36) + in(16) • (n A dtu) (46) 
+ |n(lc) ■ (n A dtii) (2c) + n(2c) ■ (n A dtn) (3c) + n(4c) ■ (n A dtii) (3c) + |n(lc) ■ (n A dtii) (4c) 
+ 



+ |n(l/i)-(nAajn)(2/i) + |n(2/i)-(nAatn)(3/i) + in(4/i)-(nAatn)(3/i)+n(l/i)-(nAatn)(4/i)} 
(A.13) 

In the above equation la, 2a, 3a and 4a denote the vertices for the subvortex a in the 
anti-clockwise sense , as we already have in the case of 1-vortex. Similarly for the other 
subvortices. It can be shown after a long calculation that Eqn. (A.13) can be written in 
the form A + B , as we had in the case of 1-vortex , where B consists of terms linear in 
712(2) or containing odd powers of 7^2(2). Thus A + B goes over to A — i? as we go from 
vortex to anti-vortex. 

In the case of a 4-votex [Fig. 4] we have a again a single lattice point as in the case of the 
2-vortex. By similar reasons the 4-vortex construction in terms of elementary plaquettes, 
breaks down in the flattened meron configuration limit. 
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